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ABSTRACT 

We present a Spitzcr IRAC and MIPS survey of NGC 2451 A and B, two open clusters in the 50-80 
Myr age range. We complement these data with extensive ground-based photometry and spectroscopy 
to identify the cluster members in the Spitzer survey field. We find only two members with 8 micron 
excesses. The incidence of excesses at 24 microns is much higher, 11 of 31 solar-like stars and 1 of 7 
early-type (A) stars. This work nearly completes the debris disk surveys with Spitzer of clusters in 
the 30-130 Myr range. This range is of interest because it is when large planetesimal collisions may 
have still been relatively common (as indicated by the one that led to the formation of the Moon 
during this period of the evolution of the Solar System) . We review the full set of surveys and find 
that there are only three possible cases out of about 250 roughly solar-mass stars where very large 
excesses suggest that such collisions have occurred recently. 

Subject headings: (Galaxy:) open clusters and associations: individual (NGC 2451); (stars:) circum- 
stellar matter; (stars:) planetary systems: protoplanetary disks; infrared: stars 



1. INTRODUCTION 

Studying the evolution of disks around young stars of- 
fers the opportunity to understand the formation of plan- 
etary systems and the process that leads to the formation 
of terrestrial planets. Protoplanetary disks accrete their 
gas content and evolve into debris disks (produced by 
the collision of larger sized bodies ) on a timescale of 10 
million years. The signatures of debris disks appear as 
mid-IR excess emission. The fraction of stars displaying 
mid-IR excess and the magnitude of the excess appears 
to decrease with time. On top of this overall trend, a 
few individual stars show extremely large excesses even 
at several hundred million years. Determining the na- 
ture and frequency of the events producing those excesses 
promises to provide an important perspective on plane- 
tary system evolution. 

Th e 30 - 150 M yr r ange is especi a lly in teresting 
since ICanupl ()2004f ) and iTouboul et all (|2007| ) showed 
that the Earth-Moon system was formed by a collision 
about 30-150 Myr a fter the Solar System was formed. 
iGorlova et al.l (j2007l ) called attention to a ~35 Myr old 
system with a large amount of debris that might have 
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resulted from s uch a giant collision [Rhee et al.l (| 2008fl: 
'Gorlo va et all (|2004l ): iRieke et al.l (|2005| ) and ISu et all 
(2006) also identify young debris systems with large 
24/im excesses. 

NGC 2451A and B are two young open clusters pro- 
jected on each other along the same line of sight that 
are within this interesting age range. Several attempts 
have been made to separate the two clusters and to de- 
termine the physic al parameters of each one of them. 
iPlatais et al.l (|2001[ ) analysed photometric and spectro- 
scopic data and used proper motions and radial veloci- 
ties to select members of NGC2451A. They fit theoret- 
ical isochrones to the cluster color-magnitude diagram 
(CMD) to calculate its distance, reddening, and age. 
They derived d= 188 pc, E(B-V)=0.01 and age=60 Myr. 
iHiinsch et al.l (j2003f ) carried out an X-ray study of the 
two clusters and identified 39 members of the A and 
49 of the B cluster using combined X-ray and optical 
data. They derived distances of 206 pc and 370 pc, 
with corresponding ages of 50-80 Myr and 50 Myr for 
NGC2451A and B, respectively. iHunsch et all |2004f ) 
completed the X-ray study with high-resolution spec- 
troscopy and refined the membership of the two clus- 
ters. The m ost recent distance and age estimate was 
performed bv lKharchenko et al.l (|2005l ). They estimated 
distances of 188 pc and 430 pc and ages of 57.5 Myr and 



2 



Balog et al. 





Fig. 1. — False color image of the center of NGC 2451. Blue: 4.5 fim; Green: 8 fim; Red: 24 fim 



75.9 Myr for NGC2451 A and B, respectively. However, 
their analysis was based on only a handful of member 
candidates (26 and 11 for NGC 245 1A and B, respec- 
tively) in co ntrast with 'Hiinsch et all (j2003f ). Therefore 
we adopt the lHiinsch et al. (2003) distances and ages and 
reddening (E(B-V) of 0.01 and 0.05 for NGC2451 A and 
B, respectively) for our analysis. 

In this paper we report on a Spitzer IRAC and MIPS 
survey of the central region of the clusters to study the 
incidence of infrared excesses in the 50-80 Myr age range. 
We supplement our Spitzer observations with a large- 
scale spectroscopic and optical photometric program to 
refine and extend our knowledge of the membership of 
these clusters. 

2. OBSERVATIONS 

2.1. Spitzer /IRAC and MIPS data 

Observations of N GC2451 were ob tained on 2004 
March 02 with IRAC (|Fazio et al.ll200l . The IRAC sur- 
vey covers about 0.6 square degree centered on HD 63032, 
a K2.5Ib star. The 12s high-dynamic-range mode was 
used to obtain two frames in each position, one with 0.4s 
exposure time and one with 10.4s. The observation of 
each field was repeated twice with a small offset, provid- 
ing 20.8s integration time for each position. The frames 
were processed using the Spitzer Science Center (SSC) 
IRAC Pipeline vl4.0, and mosaics were created from the 
basic calibrat ed data (BCD) frames u sing a custom IDL 
program (see iGutermuth et al.l (|2008f ) for details). Due 



to the 7 arcmin offset between channel 1/3 and channel 
2/4 the total area covered in all four channels is about 
0.45 sq degree. 

Source finding and aperture photometry on these im- 
ages were carried out using Phot Vis version 1.10, which 
is an IDL-GUI based ph otometry visualization tool (see 
IGutermuth et al.l (|2004[ ) for further details on Phot Vis). 
The radii of the source aperture, and of the inner and 
outer boundaries of the sky annulus were 2.4, 2.4 and 
7.2 arc-second, respectively. The calibration used large 
aperture measurements of standard stars. The zero point 
magnitudes of the calibration were 19.6642, 18.9276, 
16.8468, and 17.3909 corresponding to zero point fluxes 
of 280.9, 179.7, 115.0,64.13 Jy for channels 1, 2, 3, and 4, 
respectively. Aperture corrections of 0.21, 0.23, 0.35 and 
0.5 mag were applied for channels 1, 2, 3, and 4 to ac- 
count for the differences between the aperture sizes used 
for the standard stars and for the NGC 2451 photometry. 

We accepted as good detections those with photomet- 
ric uncertainties less than 0.2 mag, which allowed limit- 
ing magnitudes of 18.2, 17.5, 15.4 and 14.8 at 3.6 /zm, 
4.5 ^m, 5.8 /^m, 8.0 ^m, respectively. We detected al- 
most 30000 sources at 3.6 ^m, more than 20000 at 4.5 
/xm, and over 4000 and 3000 in the 5.8 fim and 8.0 /im 
images. There are more than 2500 sources that were 
detected in all 4 channels. Based on the 2MASS mag- 
nitudes of the detected sources, the limiting mass of our 
survey is < 0.1 Mq for NGC2451 A and around 0.15 Mq 
for NGC2451 B. 
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The MIPS ijRieke et al.ll2004[ ) 24 /Ltm survey was ob- 
tained on 2004 Apr. 11. It covers about 0.8 square de- 
grees centered on the same position as the IRAC survey. 
The observations were carried out in scan map mode. 
The fram es were proces sed using the MIPS Data Anal- 
ysis Tool ([Gordon et al.. .2005,) . We used Phot Vis to ex- 
tract the sources in the 24 /xm image and PSF fitting in 
the IRAF/DAOPHOT package was used to obtain pho- 
tometry. We used a 0.59 mag aperture correction and 
a 7.17 Jy zero point to convert our PSF photometry to 
the 24 fj,m magnitude scale. Similarly to the IRAC pho- 
tometry we also discarded sources with errors larger than 
0.2 mag, leaving 2625 sources in our MIPS sample. The 
limiting magnitude of our 24 /xm survey is 11.8, which 
corresponds to masses of 0.55 Mq (early M spectral type) 
and 0.9 Mq (early K spectral type) for stars with pure 
photospheric emission at the distance of NGC2451 A and 
B, respectively. 

Fig [T] shows the false color image of the center of 
NGC 2451. 

2.2. Optical photometry 

The optical observations were made with the SITe 
2048-#6 CCD camera attached to the 1.5 m telescope 
at CTIO on 2003 Jan. 22, 23 and 24. The camera 
was mounted at the f/13.5 focal position, covering a 
15 X 15 arcmin^ field-of-view with a resolution of 0.43 arc- 
sec/pixel for the entire 2048 x 2048 pixcP area. The ob- 
servations were made through Johnson-Cousins UBVRI 
filters, from the Tek #1^ filter set. 

The whole cluster was covered by 4 x 3 — 12 CCD 
frames centered on and around the brightest inner area 
at R.A. = 07 : 45 : 20, DEC = -38 : 02 : 00. One 
off-cluster area (separated by ~ 1.5 deg from the clus- 
ter center) was also imaged to sample the galactic fore- 
ground/background object population in the same line- 
of-sight. Each field was imaged three times through the 
same filter. One frame was obtained with a short expo- 
sure time (10 s for U and 5 s for BVRI) and the other 
two frames were taken with a longer one (250 s for [/, 70 
s for B and 50 s for VRI). 

The reductions of the raw frames were performed with 
standard routines using IRAF^ . After trimming the edges 
of the frames and subtracting the bias level from each 
image, the frames were divided by a master flat field 
image obtained by median combining the available flat 
field frames for each filter. Both dome flats and sky flats 
were taken at the beginning of each night and combined 
together into the master flat frames. After flat fleld divi- 
sion, the two long-exposure frames corresponding to the 
same filter were averaged to increase the signal-to-noise. 

The photometry of the cluster frames was computed 
via PSF-fitting using DAOPHOT implemented in IRAF. 
A 2nd order spatially variable PSF (varorder=2) was 
built for each frame to compensate for distortions of 
the PSFs due to either the optical imaging artifacts 
in the large field-of-view, or guiding errors that oc- 
curcd randomly on a few frames. For the PSF-model, 
f unction=auto was selected, i.e. each built-in PSF- 
function was computed and the best-fitting one was 

^ http:/ /www. ctio.noao.edu/instruments/filters 
^ IRAF is distributed by NOAO which is operated by the As- 
sociation of Universities for Research in Astronomy (AURA) Inc. 
under cooperative agreement with the National Science Foundation 



adopted automatically. In most cases the penny 1 or the 
penny2 function was the best-fitting one. The PSF-stars 
were selected interactively from a sample of the ~ 50 
brightest, non-saturated stars on each frame, omitting 
the ones with suspicious profiles and/or detectable neigh- 
bors within r = 15 pixels. The f itrad parameter was set 
as 5 pixels corresponding roughly to the FWHM of the 
frames with the lowest quality. The detection threshold 
was fixed at the 4(t level on each frame. 

Inspecting the results of the PSF photometry revealed 
that there is a small, but significant, difference between 
the results of the PSF and aperture photometry of the 
same stellar field. This difference was about the same 
(within ^ 0.03 mag) for the bright stars in one field, but 
its amount varied from field to field, ranging from a few 
hundredths up to ^ 0.2 mag. The aperture photometry 
was computed with Tap = 8 pixels aperture radius (about 
3.4"). The local sky level was estimated as the mode 
of the pixel distribution within an annulus having inner 
and outer radii of 10 and 20 pixels, respectively, cen- 
tered on each object. Because rap was more than twice 
the stellar FWHM for most frames, much of this differ- 
ence was probably due to the distorted stellar PSFs men- 
tioned above, plus the problems in determining the local 
sky level during the PSF-fitting. Therefore, we fixed the 
zero point of the photometry by selecting a bright, un- 
saturated, uncrowded reference star on each frame, and 
using its magnitude from aperture photometry as a refer- 
ence level. Thus, the final instrumental magnitudes were 
determined as follows: first, differential magnitudes were 
computed between the PSF-magnitudes of all stars and 
the reference object, then the magnitude of the reference 
object from aperture photometry was added to each dif- 
ferential magnitude. For the short-exposure frames, the 
same procedure was applied, but for the reference ob- 
ject the magnitudes from the long-exposure frames were 
selected, thus ensuring a common zero point for the pho- 
tometry of the long- and short-exposure frames. 

The transformation of the CTIO instrumental magni- 
tudes into the standard Johnson-Cousins system was per- 
formed via observat i ons o f Landolt photometric standard 
sequences (|Landoltl[l992l ). The L98, LlOl and RU152 
standard fields were imaged through each filter, weather 
and time permitting, on each night at different airmasses 
(up to 4 different elevations in the case of the L98 field). 
Since photometric accuracy is essential for the conclu- 
sions of this paper, we give the details of our calibration 
and standard transformation below. 

The following formulae were adopted for converting the 
instrumental magnitudes into standard ones: 

V-{v~ KX) = aci ■ CI + Cv, (1) 

CI = Pci ■ {ci - k,,X) + Cci, (2) 

where v and ci are the instrumental V^-magnitude and 
color index (u — b, b — v, v — r or v — i), V and CI are 
the standard magnitudes and colors, X the airmass, the 
fcs are the extinction coefficients, aci and f3ci are the 
color-dependent transformation slopes (the color terms), 
and the ^s are the zero points. 

Because the first night was nearly photometric, we de- 
cided to determine the color terms and the extinction 
coefficients from the magnitudes of standard stars ob- 
served that night. The instrumental magnitudes of these 
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TABLE 1 

Transformation coefficients derived from Landolt 
standards observed on 2003 jan. 22. 



Coefficient 



V 



(U-B) (B-V) (V-R) (V-I) 



k 
a 

P 

c 



0.15 



-2.457 



0.15 

1.064 
-1.848 



0.13 
0.035 
0.880 
-0.133 



0.00 0.05 
0.031 

1.055 1.021 
-0.083 0.822 



TABLE 2 

Transformation zero points for the three nights 



Date 


Cv 




Cb-v 


Cv'-H 


Cv-i 


2003 Jan 22 
2003 Jan 23 
2003 Jan 24 


-2.457 
-2.470 
-2.486 


-1.848 
-1.770 
-1.768 


-0.133 
-0.133 
-0.134 


-0.083 
-0.083 
-0.082 


0.822 
0.800 
0.802 




8600 
wavelength (A) 



Fig. 3. — Observed stellar spec tra (light blue/gray lines) and 
the best-fit synthetic data from the lMunari et aLl II2005I) spectrum 
library (black lines). 



stars were computed from aperture photometry using the 
same aperture and annulus radii as mentioned above for 
the cluster frames. The transformation coefficients were 
then derived by fitting Eq.l and 2 to the instrumental 
and standard magnitudes via ^^-minimization. Table [1] 
lists the parameters for each filter combination (note that 
both {B — V) and {V — I) were applied as the color index 
in Eq.l). 

With the parameters listed in Table [1] determined, the 
color terms and extinction coefficients were kept fixed 
and only the zero points were fitted to the data for stan- 
dard stars observed on the next two nights. Because 
we experienced slightly variable transparency during the 
second night, we have not attempted to adjust the extinc- 
tion correction coefficients to achieve the best fit. The 
variation of the zero points (collected together in Table 
121) is under 0.01 mag except for {U — B), illustrating the 
stability of the photometric conditions at the CTIO site 
during our observing run. 

Inserting the coefficients from Table [T] and O into Eq. 1 
and 2, the photometric data for the stars in the cluster 
fields have been transformed into the standard system. 
This resulted in a final sample of ~ 29, 000 stars having 
calibrated photometry in at least the 

V and / filters with uncertainties less than G.f magni- 
tude. 

The quality and the stability of the whole photometry 
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Fig. 2. — Spectroscopic sample on the V vs V-K diagram. Large 
dots; targets with spectra, small dots: all stars with available V 
and K magnitudes in the whole cluster area 



including the standard transformation has been verified 
by compari ng our standa r d V a nd B-V magnitudes with 
those from iPlatais et"al1 (|2001h for the 1800 stars in 
common between the two samples. The agreement is very 
good. The mean difference between the two samples is 
0.001 ± 0.06 mag in V and -0.007 ± 0.04 mag in B-V 
(the errors are the RMS of the residuals), thus, statisti- 
cally insignificant. Therefore our new CCD photometry 
in the field of NGC 2451 results in a sample containing 
29, 000 stars with calibrated data. This nearly doubles 
the amount of photometric data available for this field. 

2.3. Radial velocity survey 

We acquired AAOmega data using the Anglo- 
Australian Telescope at Siding Spring, Australia on three 
nights, February 15, 16 and 17, 2008, in acceptable condi- 
tions (clear skies with 1.5-2.5 arcsec seeing). In the blue 
arm we used the 2500V grating, providing A/ A A — 8000 
spectra between 4800 A and 5150 A. In the red arm 
we used the 1700D grating that has been optimized for 
recording the Ca II infrared triplet region. The red spec- 
tra range from 8350 A to 8790 A, with A/AA = 10000. 
This setup has the highest spectral resolution available 
with AAOmega, suitable to measure stellar radial ve- 
locities. In total, we acquired eleven field configura- 
tions centered on the open clusters. We selected the 
targets for our spectroscopic campaign based on their po- 
sitions in the optical-near-IR color-magnitude diagram: 
we targeted each source that is brighter than the 80 Myr 
isochrone plus 0.2 mag placed at the distance of the more 
distant cluster in the V vs. V-K diagram. We show the 
spectroscopic candidates in Fig. [Altogether with all the 
stars with available V and K photometry in the field. Our 
spectroscopic coverage is not complete for two reasons: 
(i) we concentrated our efforts to the central part of the 
cluster where we had Spitzer data; (ii) there were faint 
objects for which the low S/N of the spectra prevented 
measuring meaningful radial velocities. 

The spectra were reduced using the standard 2dF data 
reduction pipeline. We performed continuum normaliza- 
tion for the stellar spectra using the IRAF task oned- 
spec. continuum and then cleaned the strongest skyline 
residuals using linear interpolation of the surrounding 
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Fig. 4. — Left panel: EW of H/3 vs radial velocity. A and B show the clusters' mean radial velocities from iHiinsch et al.l | |2004| ). Left 
panel: all stars, right panel: radial velocity histogram of the two clusters using stars that show evidence of chromospheric activity 



continuum. 

Atmospheric parameters and radial velocity were de- 
termined for each star in an iterative process, which 
combined finding b est-fit synthetic spectrum from the 
iMunari et al.l ()2005[ ) spectrum library, with fitting, 
and cross-correlating the best-fit model with the ob- 
served spectrum to calculate the radial velocity. This 
approach is very similar to that adopte d by the Radial 
Velocity Experiment (RA VE) project (jSteinmetz et al.1 
120061 : IZwitter et al.l |2008[ ) , and our analysis is based on 
the same synthetic library as RAVE. Earlier results for 
other star clusters wi th the same instrument and anal- 
ysis can be found in iKiss et all (|2007l . 120081 ). Because 
of the wide range of temperatures (and hence spectral 
features), we needed three subsequent iterations to con- 
verge to a stable set of temperatures, surface gravities, 
metallicities and radial velocities. The latter are believed 
to be accurate within ±1-2 km for the cooler stars 
and ±5 km s~^ for the hotter stars in the sample (the 
boundary is roughly at 8000-9000 K). These values have 
been estimated from Gaussian fits of the cross-correlation 
profile using the IRAF task rv.fxcor and should only be 
considered as representative numbers. The atmospheric 
parameters {T^g, [AI/H], log g), due to their degeneracy, 
are accurate to about ±300 K and ±0.5 dex, with strong 
correlations (see IZwitter et al . (2008) for a thorough dis- 
cussion of the RAVE error analysis). 

To illustrate the difficulties one faces when analysing 
cool and hot stars together in the Ca II triplet region, 
we show sample spectral fits in Fig. |3l Since the Ca II 
lines almost exactly coincide with hydrogen lines in the 
Paschen series, we found that it was absolutely crucial 
to have the best-match template for cross-correlation. A 
slight template mismatch can easily lead to radial veloc- 
ity shifts of several km s~^ at this intermediate spectral 
resolution and hence one has to be very careful to opti- 
mize template selection (contrary to the commonly used 
practice in the optical range that the same template is 
used across a range of spectral subtypes or even types) . 
It is also unavoidable that as soon as the temperature 
reaches about 9000 K, the broad spectral features will 
lead to a degraded velocity precision simply because of 
the broadened cross-correlation profile. NGC 2451A and 
B, as reasonably young open clusters, still host a sig- 
nificant number of hotter main sequence stars, possibly 
resulting in degraded velocity precision for a fraction of 
stars. 



Another difficulty was the very high incidence of Ca II 
emission (defined by the excess flux of the calcium lines 
relative to the best-fit model spectra), especially in late- 
type main sequence stars in the clusters. Both clusters 
are young, hence the high rotation rates can lead to el- 
evated chromospheri c activity, for which t he Ca II lines 
are good indicators (jAndretta et al.|[2005f ). When com- 
paring the best-fit models to the observed spectra, we no- 
ticed distorted Ca II profiles for many stars, frequently 
in the form of central emission components within the 
lines. These can introduce random velocity errors of up 
to 10-20 km s^^ in the maximum of the cross-correlation 
profile. This was particularly apparent after modifying 
the velocity determination for 405 stars with Ca II line 
profile irregularities: after excluding the Ca II lines from 
cross-correlations, a large fraction of the refined velocities 
were very close to the cluster's mean velocities (FigH]). 

We also measured cross-correlation velocities from the 
lower resolution blue spectra. These data were less useful 
because many of the stars have a broad II/3 line and only 
a few additional weak features in the recorded wavelength 
range, leading to a broad cross-correlation profile and 
radial velocities accurate to only about ±10 km s~^. 

3. ANALYSIS 

3.1. Radial velocities of the clusters 

The clusters are in the plane of the Milky Way so the 
background contamination overwhelms the cluster mem- 
bers. To estimate the radial velocities of the clusters we 
need to make some assumptions to separate the proba- 
ble cluster members from the field based on our present 
knowledge of the cluster. Based on the young age of 
the clusters, we assumed that our sample should contain 
a large number of K and M dwarfs that show chromo- 
spheric activity, i.e., signs of emission in the CaT lines. 
We identified 405 active red objects in our sample. Sev- 
eral of our target stars show sharp H/3 emission as addi- 
tional evidence for activity. In the left panel of Fig. [4] we 
show the equivalent width of H/? vs radial velocity (nega- 
tive EWs mean emission). We see two distinct groups of 
emission line objects in the figures around the approxi- 
mate position of the radial velocity values of NGC2451 A 
an d B (22.7 kms~ ^ and 14 kms~^ ) that were estimated 
by IHiinsch et al.l ([2004) . We show the radial velocity 
histogram of these stars in the right panel of Fig. [H 
From this histogram we calculate mean radial velocities 
26.7±2.4 kms-i and 17.4±2 kms^^ for NGC2451 A and 
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20 40 60 

Fig. 5. — Radial velocity histogram of the two clusters using all 
stars with log g > 3.0. A and B show the radial velocity values de- 
termined using the chromospherically active young stars for NGC 
2451 A and B, respectively. 

B, respectively. Because chromospheric activity is not 
always present in young stars (e.g. early type ones) we 
use this test only to estimate true cluster velocities. 

Our average radial veloci ty value s are a bout 3 km s~^ 
higher than those of Hiins chet al.l (l2004f). Considering 
the formal standard deviation of the' Hiinsch et al.l ()2004l ) 
data (2.3 km s~^ for cluster A and 2.5 km s~^ for cluster 
B) and the small number of member stars in that sample 
(21 for cluster A and 10 for cluster B), the 3 km s~^ 
systematic difference is hardly significant. Given that 
our sample contains three-to-ten times more members 
for each cluster, we adopt our mean values for the cluster 
velocities. 

3.2. Cluster membership and color magnitude diagrams 

In Fig. Owe show the radial velocity histogram of our 
complete sample with log g > 3 (giants filtered out) . The 
two peaks at the positions of the radial velocities of the 
two clusters (designated with A and B) are readily dis- 
tinguishable. However a large number of contaminating 
stars might be present due to the strong galactic back- 
ground so the radial velocity is not enough to determine 
accurate membership status. Therefore, we used photo- 
metric diagrams also to clean our initial sample. 

We performed the final determination of memberships 
as follows. We selected member candidates using the 
radial velocity and log g criteria. We accepted all objects 
as cluster members that had log g > 3 (1079 out of 2570 
objects satisfied this criterion) and not more than a 2 cr 
deviation from the clusters' mean radial velocities (13.4 
kms-i < Yrad < 21.4 kms'^ for NGC2451 B and 21.9 
kms-i < Yrad < 31.5 kms"! for NGC2451 A, 152 and 
200 objects for NGC2451 A and B, respectively). Then 
we plotted the initial sample in color magnitude diagrams 
and removed the object s with posit i ons in consistent with 
the 50 Myr isochrone of lSiess et al.l (|2000f ) at the distance 
of the clusters. We show the final cleaned diagrams in Fig 
[6l Our final sample consists of 60 objects in NGC2451 A 
and 121 objects in NGC2451 B The photometry of the 
members is listed in Tables [Sj SI [SJ and [H 

3.2.1. Members from \Platais et al.\ (200A) 

There are 18 me mbers of NGC2451 A reported by 
iPlatais et al.l (j2001) in the area covered by all 5 bands of 
our 5piteer observations. Thirteen of them were included 
in our radial velocity survey. We rejected 7 of these stars 



based on their radial velocity; 4 were included in our 
clean NGC2451 A membership list while 2 have radial 
velocities just outside the 2 a limit. We included these 
two stars along with the other 5 not covered by the ra- 
dial velocity survey in our final sample for NGC2451 A 
(they are marked with filled square s in the plots in the 
left column of FigO. We used the IPlatais et al.l (|200lD 
magnitudes for these objects because they were not de- 
tected by our UBVRI survey. As a result, they are not 
plotted in the V vs V-I plot in FigjSj 



3.2.2. Members from\Hunsch et aV [2003) and 
\Hunsch et al. (2004) 

iHiinsch et all (|2003[ ) hst 87 objects in the field of 
our Spitzer survey as candidate members of either 
N GC2451 A or B . Unf ortunately the coordinates given 
m iHunsch et all (|2003D are from X-ray data and have 
sufficiently large uncertainties that it is not easy to iden- 
tify their sources in the IRAC and MIPS datasets. We 
could identify 58 sources out of 87 with a detection in 
at least one Spitzer band. Thirty four of them were in- 
cluded in our radial velocity survey which also inc l uded 7 
stars identified as non-members by |Htinsch^^2j[200^ 

In the case of NGC 2451 A, sixteen iHiinsch et al.l 
(|2003f) members are included in our radial velocity sur- 
vey, which confirms t he membership statu s of 10. In 
NGC2451 B we have 13 IHiinsch erall ()2003D members in 
our spectroscopic sample. Eight of them are confirmed as 
members by the radial velocity data while two have radial 
velocities indicating th at they are memb ers of NGC2451 

A. We can reject 3 Hii nsch et al.l ()2003l ) members based 
on their radial velocit y data. There is a third group 
of IHiinsch et al.l (|2003f l stars that contains objects that 
might be members of either of the clusters. Six of them 
are included in our radial velocity survey. Based on the 
radial velocity data we confirm membership of 5 objects 
(4 in NGC2451 A and 1 in NGC2451 B). We included 
the iHiinsch et al.l ()2003l ) members without radial veloci- 
ties in our final sample if they were detected in at least 
one of the Spitzer bands (3 and 10 for NGC2451 A and 

B, respectively). They are marked with filled triangles in 
Fig [6l These additions make the total number of cluster 
members 69 and 131 for NGC2451 A and B, respectively. 
We detect 22 members at 24/im, 47 at 8 and 4.5/im, and 
49 at 5.8 and 3.6 fim in NGC2451 A. In NGC2451 B we 
detect 20 members at 24^m, 92 at 8 //m, 96 at 5.8 ^m, 
94 at 4.5 fi, and 98 at 3.6 ^m. Our detection rate in 
the IRAC bands is about 70% in both clusters while at 
24 /im it is 32% in NGC2451 A and 15% in NGC2451 
B, resp ectively. These detec tion rates are very simlar to 
that of iGorlova et all ()2007l ) for the --35 Myr old cluster 
NGC2457. In the following we use only the sources that 
are detected in all four IRAC bands or detected at 24/im. 

4. IDENTIFYING EXCESS CANDIDATES 

4.1. IRAC excesses 

The four IRAC bands covering the 3.6—8 /im regime of 
the spectrum trace dust from the near-sublimation tem- 
perature of 1000 to 400 K, corresponding to radii within 
1 AU for solar- type stars. We can identify 8 /im excesses 
from the IRAC data, if we plot K-[8] vs V-K. Fig|7|shows 
the V-K vs K-[8] color-color diagram of NGC2451 A (left 
panel) and B (right panel). The large baseline of the V- 
K color makes this filter combination an ideal choice for 
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Fig. 6. — Color-magnitude diagram of the members of the two clusters (left panel: NGC2451 A; right panel: NGC2451 B). dots: stars in 
our spectroscopic su rvey, squares: members from Plat ais et al, {2001) not included in our spectroscopic sample, triangl es: members from 
IHiinsch et al.l i2003l) not included in our spectroscopic sa mple. The solid line show the 50Myr and 80Myr isochrone of ISiess et al. | l|2000|) 
for NG C2451 A and the 50 Myr isochrone of Siess et all II2000I ') for NGC2451 B, shifted with the distance and extinction of IHiinsch et al.l 
(|2003) 
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photometric spectral type determination given the low 
reddening of the clusters. There is a monotonic relation- 
ship between increasing V-K and stellar spectral type 
that results in a narrow locus, extending from early to 
late type stars; the scatter in this locus increases for V- 
K>5 due to measurement errors. Objects to the right of 
this locus are stars with excesses at 8 fim. 

The robust determination of cluster membership allows 
us to establish a good color-color relation for 50-80 Myr 
clusters by fitting a 2nd degree polynomial to the locus of 
the V-K - K-[8] and V-K - K-[5.8] diagrams. We merged 
the data of the two clusters to increase the number of 
sources in the fit and required both colors of an AO star 
to be 0. Eq. 3 and 4 show the results of the fit which is 
also shown with a solid line on Fig [7] and [H The rms of 
the datapoints around the fit is 0.05 in both colors. 

(Ks - [8])o = 0.0115 ± 0.0011(V^ - Ks)l 

-5.89 X 10"^ ± 0.0053(^^ - Ks)o (3) 
{Ks - [5.8])o = 0.0125 ± 0.0012(1/ - Ks)l 

-0.0078 ± 0.0060(1/ - Ks)o (4) 

No star with 8 /im excess (relative to K) was detected 
in NGC2451 A and only two candidate excess stars (#1 
and =ff2) and another two weak possible excess candidates 
(#3 and #4) were detected in NGC 2451 B. We also plot 
the same diagrams but using the IRAC 5.8 /im channel 



on FiglHl This diagram allows us to check whether our 
8 ^m detections are real or are just a result of some 
anomaly in the image (e.g. a cosmic ray hit). Both 8 /xm 
excess candidates are confirmed by having excess at 5.8 
/im also. However star #2 has a close companion that 
might compromise the 8 //m and 5.8 /im photometry (see 
Fig [T0|) . Stars #3 and #4 are also situated slightly to 
the right of the main locus of sources in V-K vs K-[5.8], 
however, #3 has a larger excess at 5.8 fim so its status 
as an excess candidate is questionable. None of the four 
8 //m excess candidates were detected at 24/im. 

We fitted a Kurucz model to the short wavelength por- 
tion of the sources' SEDs. Fig. [TT] shows the results of 
the fit. Based on their SEDs wc can conclude that only 
one of the four candidates (#1) is a true excess source. 
The SED of #2 shows that the photometry of this source 
is contaminated by the nearby bright star because the ex- 
cess extends over all IRAC bands. The SED of #3 shows 
no excess at all and that of #4 has a marginal one. The 
significance of the excess for #4 calculated from the pho- 
tometric errors and the rms of the fit is slightly lower than 
3cr. 

We noticed that there is a small hint of an 8 fim excess 
on the SED of the 24/nm excess star #12 (see next sec- 
tion) and all but one of the A (or earlier) type stars have 
K-[8] colors less than zero although they are inside the 
3(7 limit of the trendline in Fig [71 To check the origin of 
this slight discrepancy we replaced K with [3.6] and used 
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the V-K vs [3.6]-[8.0] diagram to detect excesses (Fig. 
[9]). A linear relation between the V-K and the [3.6]-[8.0] 
colors can be found for the photospheric sources 

([3.6] - [8])o = 0.01842 ± 0.001277(y - Ks)q (5) 

with rms scatter of 3%. We also show this relation along 
with the 3(7 limits in Fig. [9] 

In this diagram there is no systematic deviation from 
the main relation and it seems that #12 has a marginal 
excess at 8 /im relative to 3.6/im. The significance of the 
excess is about 4(t. There are some excess candidates 
among the late M type stars of NGC2451 A but their 
significance is less than 3 a. In the case of NGC2451 B 
(Fig. [9] right panel) the diagram clearly shows that the 
only viable 8 fj,m excess candidate is #1. Although #3 
shows some [3.6]- [8.0] excess, the SED of the source does 
not support it. 




Fig. 10. — Location of star #2. The green circle shows the 
aperture size that was used for the photometry 



To summarize, at 8 fim we have less than 2% excess 
fraction in NGC2451 A (one marginal excess) and a less 
than 1% excess fraction in NGC2451 B (one clean ex- 
cess) . This result is very similar to that of lGorlova et al.l 
([2007) suggesting that IR excess shortward of 24 ^m is 
very rare in clusters older than 20-30 Myr and indicates 
that the disk within 10 AU is cleared of dust for nearly 
all stellar systems older than 30 Myr. 



4.2. MIPS excesses 

Fig [H shows the dereddened V-K vs K-[24] for 
NGC2451 A (left panel) and NGC2451 B (right panel). 
Because there are too few detections to determine a zero 
excess locus purely from these data, for this purpose we 
have used a database of more than 1500 stars derived 
from nearly all of the Spitzer debris disk programs. The 
sample has good coverage from -0.4 < V-K < 3.5, with 
still a number of members up to V-K = 8. A simple poly- 
nomial fit has been made to the photospheric color-color 
locus. The rms scatter around this fit (clipping the pos- 
itive outliers that have excesses) is only 3%. The locus 
itself is therefore very well determined and we can take 
a typical uncertainty for the measurement of any given 
star to be 3% relative to this locus (plus any statistical 
errors). This locus with 3-a errors (9%) is plotted in Fig 

m 

It is obvious that excess at 24 /im is far more common 
than at shorter wavelengths. We find 9 stars in the case 
of NGC2451A and 10 stars in the case of NGC2451 B 
to the right of the photospheric locus. However, only 
7 of them in both diagrams (the numbered sources in 
Fig [T^ are more than 3a distance from the locus (the 
distance of the negative outlier left of the locus). The 
measured (K-[24])o color along with the predicted photo- 
spheric K-[24] color of these objects and the significance 
of the excesses are presented in Table [71 For the predic- 
tion of the photospheric colors we used the photospheric 
locus shown in Fig [T^ The significance was calculated 
using the photometric erors and the rms around this lo- 
cus. To check the validity of our excess candidates (as in 
the previous section) we fit Kurucz model atmospheres 
to the short wavelength part of their SEDs. The fits are 
shown in Fig [13] and [TH In all cases the 24 /im flux 
is above the photosphere by more than 3 a, confirming 
the excess. The fraction of stars with excess at 24 /im 
is about 33±10 % (7/22) for NGC2451 A and 36±10 
% (7/20) for NGC2451 B. We calculated the disk frac- 
tions and errors usin g Bayesian statistics as described in 
iGaspar et all (|2009D . All but one of our 24 fim excess 
candidates have spectral type F or later (based on their 
V-K colors and fitted SEDs). We detect only one star 
earlier than FO with 24 /^m excess in NGC2451 A and 
none in NGC2451 B. 

Our excess fractions among solar type stars are 35±12 
% (5/15) in NGC251 A and 39±12 % (6/16) for NGC 
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2451 B, and for early type stars are 29±17% (1/5) and 
2G±20 % (0/3) for NGC2451 A and B, respectively. 



5. DISCUSSION 

5.1. Excess fractions 

Our results for excess fraction in the solar type stars 
of NGC2451 B are roughly consistent with other clus- 
ters with sim i lar ag es such as IC 2391 (50 Myr - 31% 
200^ ;^ or NGC 2547 (35 Myr - 40%; 



Siegler et al.l 



Gorlova et al.l 120071) ). In the case of NGC2451 A our 
excess fraction is significantly large r than that of M47 
(80 Myr - 6% (jGorlova et al.ll2004f) ) which might indi- 
cate that the age of NGC2451 A is closer to 50 Myr 
rather than 80 Myr. However, the study of M47 has a 



higher threshold for detection of excesses, which might 
also account for the lower fraction. 

For the earl y type stars o ur res ults are in agr e ement 
with those of ISiegler et"aLl (|2007D and ISu et al.l (|2006D 
who found a very low excess fraction ( about 10% and 
14% for IC2391 (50 Myr) and IC2602 (30 Myr), respec- 
tively) . 

Fig [T5] shows the d isk fractions a s a fu nction of age 
(similar to Fig 9. of iGaspar et al.l (|2009D ) (left panel: 
early type stars; right panel: late type sta rs). Open cir- 
cles sh ow the data from Table 4 and 6 of iGaspar et al.l 
([20091 ) while filled squares designate the two new data- 
points from our study. We note that for stars younger 
than 10 Myr the excess not necessarily indicates the pres- 
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Wavelength [microns] 



Fig. 13. — The SEDs of the 24/im excess candidates in NGC2451 A together with the fitted Kurucz atmospheric model. The errorbars 
represent the 3 a photometric error. 



ence of debris but primordial dust lef t over from t he sta r the debris disk fraction. An interesting feature of the fig- 
formation process as pointed out by iRhee et"al] (|2007f ). ures is that there is a dip in the fraction of debris disks in 
so the first three points in the figure might overestimate the age range of 30-80 Myr. Outside this range the decay 
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Fig. 14. — The SEDs of the 24;um excess candidates in NGC2451 B together with the fitted Kurucz atmospheric model. The errorbars 
represent the 3 a photometric error. 



seems to be monotonic. Our two clusters slightly deviate main locus of stars along with some other clusters of sim- 

from the main trend in the late type stars also but in ilar age. Incompleteness might account for the discrcp- 
a different direction. They are both situated above the ancy. However, even if we remove the faintest stars (#17 
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TABLE 7 

Measured and predicted photospheric (K-[24])o colors of cluster members 

WITH 24 excesses 



ID 


RA 


DEC 


(K-[24])o 


<7{(K-[24])o 


) iK-mf' 


significance of 




J2000 


J2000 


mag 


mag 


mag 


excess 




NGC2451 A 


#6 


116.04942 


-37.73994 


1.099 


0.084 


0.271 


9(7 


#7 


116.74747 


-38.37820 


0.648 


0.063 


0.023 


9cr 


#8 


116.31223 


-38.19624 


0.547 


0.077 


-0.018 


la 


#9 


116.23017 


-38.07956 


0.574 


0.041 


-0.019 


Via 


#10 


116.53718 


-37.39910 


0.268 


0.043 


-0.018 


ha 


#11 


116.10025 


-38.59977 


0.235 


0.042 


-0.011 


5o- 


#12 


116.18701 


-38.05376 


1.604 


0.040 


-0.024 


33cr 


NGC2451 B 


#13 


116.67488 


-37.98390 


1.422 


0.076 


0.353 


13cr 


#14 


115.95562 


-37.89574 


0.799 


0.040 


-0.022 


Via 


#15 


116.21267 


-37.66869 


0.496 


0.041 


-0.023 


\\a 


#16 


116.50776 


-38.00752 


0.516 


0.047 


-0.022 


lOcr 


#17 


116.44107 


-37.89720 


0.493 


0.052 


-0.010 


9a 


#18 


116.74249 


-37.63916 


0.293 


0.061 


-0.024 


5o- 


#19 


116.36051 


-38.19643 


0.317 


0.043 


-0.025 


la 
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Fig. 15. — The fractio n of star with 24/xm excess vs age in open clusters. Left panel: early type stars. Right panel: solar type stars. 
Open circles: data from lIGaspar et al.ll200S) : filled squares: NGC2451 A and B 



in Figini) in NGC2451 B and the three faintest ones (#7 
and the two photospheric sources around V-K~3 Fig [12] 
left panel) in NGC2451 A, the excess fractions (35±12 
% and 36±13% for NGC2451 A and B, respectively ) are 
still slightly above the main locus. 

Unfortunately low number statisics might reduce the 
reliability of these results. To check whether they are 
statistically signific ant we binned the cluster data of 
iGaspar et al.l ((2009') supplemented with our two clusters 
into five equal logarithmic age bins 3-lOMyr, 10-31. 6Myr, 
31.6-lOOMyr, 100-316Myr, 316-1000 Myr. We show the 
result of the binning in Fig [16] (left panel: early type 
stars;right panel: late type stars). There is a hint of 
a drop (about 1 ct) in the excess fraction of early type 
stars in the 31.6-100 Myr range then, a rise and another 
drop after 316 Myr. For the late type stars a small bump 
breaks the continous decay with about the same signifi- 
cance as the dip in the early type sample. More detailed 
and very accurate membership study of the clusters in 
this age range is necessary to decide whether this behav- 
ior can be attributed solely to statistical fluctuation. 



5.2. Stars with large excess 

Among all stars in NGC 2451 A/B showing excess at 
24/im, #12 (HD 62938) has by far the largest; K-[24] 
= 1.6 corresponding to an excess ratio (R) more than 4 
meaning that the 24/xm flux density is more the 4 times 
larger than the expected photospheric flux level. The 
MIPS data place a 2a upper limit of 63 mJy on the flux 
density at 70 /xm, requiring a color temperature between 
24 and 70 /Ltm > 95 K. To put this star into context with 
other results we performed a literature search for main- 
sequence stars older than 20 Myr with similarly large 
or larger excess. We found 9 stars with excess ratios 
larger than 4 at 24/im. We also included two famous 
somewhat younger large excess stars HD39060 (/3 Pic) 
and HD181296 {rj Tel) for comparison. They are listed 
in Table [5] The excess of HD21362, a B6 star, is due to 
free-free emission so we rejected this star from our final 
sample. 

Fig [17] shows the excess vs age diagram f or our sam- 
le of 10 stars together with the sample of iRieke et al.l 
20051 ). We express the e xcess in 2.5xlog{R) so we can 
compare the lRieke et al.l ()2005| ) data with other publica- 
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Fig. 16. — The binned excess fraction of open cluster vs age. Left panel: early type stars. Right panel: solar type stars. 



TABLE 8 
Stars with large excess 



Name 


SpT 


Age 


excess ref. 


HD21362t 


B6Vn 


80 Myr 


Rieke et al. (2005) 


HD39060 


A5V 


12 Myr 


Rieke et al. (2005) 


HD109573 


AOV 


10-20 Myr 


Rieke et al. (2005') 


NSV17775t 


AlV 


25 Myr 


Young et al. (2004) 


NGC2547/id8 


G7V 


35 Myr 


Gorlova et al. (2007) 


NGC2547/id7 


M 


35 Myr 


Gorlova et al. (2007) 


P1121 


F9V 


80 Myr 


Gorlova et al. (2004) 


BD+20 307 


GO 


2 Gyr 


Song et al. (2005) 


HD23514 


F6 


120 Myr 


Rhee et al. ("2008) 


HD181296 


AOV 


12Myr 


Su et al. (2o3ST 


HD62938 (#12) 


AOV 


50-80 Myr 


this work 



age ref 



Sone et al. (2001 ) 
^rtega et al. (2002, 2004) 
^ckerrnan & Song (2004a) 
Jeffries & Oliveira (2005) 
Jeffries & Oliveira (2005) 
Jeffries & Oliveira (2005) 
Roio Arellano et al. (1997) 
Song et al. (2005) 
Rhee et al. (2008) 
Zuckerman fc Song (2004b) 
Platais et al. 120011 



t excess comes from free-free emission^ incorrectly listed as a 
Cataclysmic Variable Star in SIMBAD 



tions where the excess is expressed in K-[24] mag. All 
but three of the stars with lar ge excesses f a ll wit hin the 
envelope of the decay trend of i Rieke et al.l (120051 ). (The 
normalization of this envelope is somewhat arbitrary so 
no distinction should be made between star s just above or 
j ust b elow it.) Two exc eptions, P11 2 1 of iGorlova et al.l 
(|2004[ ) and HD23514 of iRhee et all (120081) are young. 
The third one BD-h20 307 (ISong et al.ll2005l) is very old 
(about 2 Gyr (jZuckerman et al.l 120081 )). The excesses 
of these stars are much higher than would be predicted 
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Fig. 17. — 24/^m excess vs a ge. Squares: s tars w ith large excess 
from Table [S] dots: stars from lRieke et aEI II2005I ). solid line: the 
normal decay trend from infrared excesses 



based on their ages. Another two stars that are older 
than 30 Myr are borderline cases close to the upper en- 
velope of the deca y trend. One of th em is a G7 star (ID8 
in NGC 2547 Go rlova et all (|2007l )) the other is an M 
star in the same cluster (1D7). We retain these two stars 
as candidate extreme excess stars. 

Although these stars were identified only on the basis 
of their large excesses at 24/im, they have a number of 
unique characteristics in comparison with other debris 
disks. For example, without full spectral energy distri- 
butions we cannot compute the fractional luminosities, 
Ld/L*, but we can put lower limits on this parameter. 
As an illustration, we assume a star represented by a 
6000K blackbody and a blackbody-like excess ten times 
the stellar output at 24 /zm. The minimum fractional 
luminosity is then 2 x 10~^ at an excess temperature of 
between 150 and 160K. Seven of the ten large excess stars 
are at this level or above. The minimum fractional lu- 
minosity for the remaining three large excess sources is 
> 2 X 10^''. For comparison, debri s disks typi cally have 
fractional luminosities below lO""*^ (jMoor et al.ii2006 ). In 
addition, many of these systems are known to have strong 
spectral features in the mid-infrared, establishing that 
their emission there is dominated by w arm, very finely 
divided dust (/3 Pic. lChen etlHI (|2007i ): NGC 2547 IDS 
and PI 121 , Gorlo va et al. in preparation ; BP + 20 307, 
ISong et all (pOOl : HD 23514. iRhee eFall (pOOl l. 

These latter sources represent extreme examples of 
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debris disks containing warm dust. Additional ex- 
amples of systenis with warm dust include Fomal haut 
(jStapelfeldt et al.ll2004D : /? Leo, (lAkeson et al. |l2009l. our 



unpubli shed observations) ; 7? Tel, (ISmith et a l. 2009); n 
Corvi, (jSmith et al .' 2008); HR 4796 A, (Wahhaj.eLaLl 



[2005h: e Eri, 
(jLisse et al.| [200 



Backman et a l. 2009); HP 



HP 69830, (Beichman et al 



and C Lev. (|Moerchen et al.ll2007( ). iMorales et al 



1 13766 , 
12001 : 



mm ) 



report 28 additional systems with dust at T > 180 K and 
that dominates the disk emission at 24 fj,. It is therefore 
plausible that many of the other large excess systems 
have substantial emission by warm dust. 

Although many of the large-excess stars are roughly 
of solar type, even in this case they are relatively un- 
common. In particular, for the 30 - 130 Myr age range. 
Table [S] shows there are only three solar-like extreme 
excess stars. They are drawn from a parent popu- 
lation of all the stars observed within this age range 
and with spectral types between mid-F (F 4) and mid- 



KJKA), a total of about 2 50 (from M47, Go rlova et al 
2004) ; the Pleiades, St auffer et all (12005). Gorlova et al 
20061 ). and our un published work ; IC 23 91. ISiegler et al 



Mm-. NGG 254 7.lGorlova et all (l200l : the FEPS sam- 
ple. lMever etlHI (12008[); and NGG 2451 A&B, this work). 



That is, the incidence of extreme excess systems among 
solar-like stars in the 30-130Myr age range is only about 
1%. 

The 30-130Myr age range insures that these extreme 
excesses are unlikely to originate in the wave of oli- 
garth plane t building that peaks around 10-15 Myr (e.g., 
iGurrie e'Ta l. (2008) and references therein). There are al- 
ternative possibilities for generating them, however. One 
is collisions between large bodies, simi lar to the collisio n 
that led to the formation of the Moon (jRhee et al.ll2008f ). 
It is thought that about three moon masses of material 
were thrown into orbit aro und th e Sun in the impact that 
formed our Moon (|Ganup|l2004D . The rapid injection of 
such a large mass would produce a signature roughly in 
agreement with the observations of these large-excess ob- 
jects (Gorlova et al. in preparation). A rough e- folding 
time of 2 million years for the dissipatio n of this mass 
can b e estimated from the calculations of iGrogan et al.l 
(|2001| ). In addition, a process that results in very vig- 
orous gravitational stirring of a dense asteroid belt (e.g., 
thr ough gravitationally resonant orbits of massive plan- 
ets, iGomes et al.l (|2005l)) is also a candidate to produce 
huge infrared excesses. 

Although Table[5]contains several early type stars older 
than 20 Myr and with large excesses, none of them is 
above the decay trend of normal debris disk stars. In 
fact, as Table |S] emphasizes, all the large-excess early 
type stars are relatively young and may just still be in 
the decay of their initial peak debris output. The place- 
ment of the trend line is rather arbitrary. Nevetheless, we 
have considered why no very large excesses are found for 
older early- type stars. We calculated what mass would 
be needed to cause an excess about 4 times above the 
stellar photosphere at 24 ^m in an A star. We used only 
grains larger than 6.3 ^m, because radiation pressure will 
eject small grains from the system, with terminal velocity 
V = ^^^iP ~ l)^ where M* is the mass of the central 
star, rinit is the radius of the initial orbit of the grain and 
P is the ratio of forces from radiation pressure and grav- 



ity. iPominik fc Peciril (|2003f ) showed that the orbiting 
grains normally dominate the visibility of the debris disk 
over the blown-out grains. To get the total mass needed 
to cause the required 24^m excess we first assumed that 
all the flux comes from the smallest grains then we cal- 
culated the total surface area of these grains. Using the 
surface area and the N{a) ^ a~^-^ grain size distribution 
we were able to determine the scaling constant using the 
following equation 



A - 

Att — ca 
4 



(6) 



where amin is the diameter of the smallest grain reman- 
ing in the system and a^ax is a nominal maximum parent 
planetesimal size (100m). Using the scaling constant we 
can calculate the total mass by integrating ^■ko?' pca~^'^ 



from a„ 



to an 



We found that the total mass needed 



to produce the excess is on the order of a couple of 
Moon masses. The same calculation for G-type stars 
with amin = 0.01/im leads to two orders of magnitude 
less mass (~ 10~^MMoon)- 

In the case of A stars, radiation pressure clears the 
small grains from the disk in a few hundred years while 
these g rains re main around the G-star for a few times 10^ 
years (Ghen e t al. 2005). To quantify the effects of this 
difference, we also calculated the scaling factor and the 
total mass neglecting the effect of blowout for the A star 
and integrating over the entire size range from 10^^//m 
to 100m. We found that in this case we would need 
more than 25 times less material for producing the same 
level of excess around the A star. In this hypothetical 
case the masses needed for A- and G-type stars are the 
same order of magnitude. This implies that their ability 
to eject grains very efficiently through photon pressure 
plays a large role in the absence of extreme excess around 
older early type stars. Those small grains carry a large 
surface area (relative to their mass) which would make 
the excess more easily detectable. A corollary is that 
large grain production rates around early-type stars will 
be indicated by large flows of small grains out of the 
system driven by photon pr essure. Possible examples 
include Vega (|Su et al.ll2005h and HR8799 (Su et al. in 
preparation) 

6. CONCLUSIONS 

We present optical and 5piteer/IRAG-MIPS photom- 
etry and medium resolution spectroscopy of the central 
0.75 sq deg. region of NGG 2451 A and B. Using ra- 
dial velocity data and optical/near-IR color magnitude 
diagrams we selected 61 members of NGG 2451 A and 
131 members of NGG2451 B. We found one object in 
NGG2451 B and none in NGG2451 A with excess at 8^m. 
We conclude that 8 micron excesses at the ages of these 
clusters are probably associated with an isolated event 
that disturbs the distribution of planetesimals to yield a 
sequence of collisional cascades producing dust. 

We find one object with 24/i,m excess among the early 
type stars of NGG 2451 B and none in NGG2451 A. For 
the solar type stars we found 5 out of 15 stars with excess 
in NGG2451 A and 7 out of 16 for NGG2451 B. In general 
these results agree with excess fractions reported in the 
literature. 

We used the data presented in iGaspar et al.l ()2009D 
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supplemented with our new datapoints to plot the ex- 
cess fraction vs age, and found that the excess fraction 
unexpectedly drops in the range of 30-80Myr for early 
type stars while it slightly rises in the same age for late 
type st ars. We tested our findings by binning the cluster 
data of iGaspar et alJ f2009,) into equal logarithmic age 
bins. The results are consistent with some structure but 
also indicate that the fluctuations may be statistical in 
nature. 

We summarize the detections to date of large debris- 
disk 24 /im excesses. There are only five extreme cases, 
with fractional luminosities > 2 x 10~^. We did not find 
extreme excesses around older early type stars, excesses 
that would indicate a catastrophic coUison between form- 
ing planetesimals. We found that blowout of the small 



grains probably plays a very important role in the lack 
of detections. 
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TABLE 3 

UBVRIJHK Photometry of the members of NGC2451 A. 
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-37.96688 
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0.009 


9.682 


0.013 


9.547 


0.023 


9.511 


0.020 


9.470 


0.019 


116.03583 


-37.98349 


9.346 


0.001 


9.466 
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0.031 
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0.020 
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116.07976 
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0.023 
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0.024 
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116.20300 


-38.09869 
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9.898 


0.022 
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0.023 
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0.004 
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0.026 
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13.43 
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11.327 


0.025 


10.983 


0.033 


10.876 


0.028 



TABLE 4 — Continued 



RA 


DEC 
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<x{U) 
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a{B) 
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a(V) 
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a(R) 
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116.39170 


-37.79194 


13.020 


0.003 


12.996 


0.005 


12.425 


0.007 


12.112 


0.004 


11.784 


0.005 


11.310 


0.019 


11.073 
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116.50776 


-38.00752 (#16) 


12.613 


0.001 


12.723 
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0.001 


11.384 


0.002 


10.863 


0.021 


10.606 


0.022 


10.515 


0.019 



numbers next to the coordinates designate 24 cxcees sources from Table [7] 



TABLE 5 

Spitzer 3.6/.iM - 24/^M PHOTOMETRY OF MEBERS OF NGC2451 A 



RA 


DEC 
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a([3.6]) 
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a([4.5]) 
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a([5.8]) 
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-38.11860 
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0.003 
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0.003 


10.363 
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0.008 






115.97250 
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0.014 
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0.023 
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-37.74447 
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13.986 
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13.812 


0.124 






116.13418 


-37.68301 


13.068 


0.008 


12.948 


0.011 


12.971 


0.044 


12.976 


0.062 






116.13475 


-38.03229 


9.382 


0.002 


9.403 


0.003 


9.374 


0.005 


9.365 


0.005 


9.245 


0.034 


116.14048 


-38.15269 
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0.007 


11.881 
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11.828 


0.035 
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-38.32344 


13.070 
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116.20787 


-38.16160 


10.252 


0.003 


10.205 
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10.222 


0.007 


10.180 


0.008 


10.161 


0.051 


116.21219 


-37.85936 


10.638 


0.003 


10.695 


0.004 


10.671 


0.009 


10.647 


0.014 


10.708 


0.067 


116.21249 


-37.74415 


13.037 


0.009 


12.987 
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12.916 


0.047 


12.941 


0.061 






116.23017 


-38.07956 (#9) 


10.033 


0.003 


9.995 


0.003 


10.038 


0.007 


9.963 


0.007 


9.465 


0.036 


116.25277 


-37.53289 






















116.25701 


-38.03892 


8.970 


0.002 


8.981 


0.003 


8.914 


0.005 


9.039 


0.005 


9.012 


0.033 


116.27029 


-38.07462 


12.226 


0.007 


12.145 


0.008 


12.139 


0.022 


12.110 


0.041 






116.28444 


-38.14332 


11.223 


0.004 


11.278 


0.005 


11.199 


0.012 


11.210 


0.018 


11.209 


0.141 


116.29981 


-37.79209 


12.534 


0.006 


12.477 


0.008 


12.413 


0.027 


12.439 


0.039 






116.31223 


-38.19624 (#8) 


9.193 


0.003 


9.162 


0.004 


9.136 


0.005 


9.137 


0.004 


8.645 


0.046 


116.33550 


-37.89964 


11.210 


0.004 


11.219 


0.005 


11.138 


0.023 


11.135 


0.014 






116.36404 


-37.59628 






7.220 


0.002 






7.189 


0.002 


7.230 


0.031 


116.38750 


-37.51380 






















116.38871 


-38.27644 


12.409 


0.006 






12.271 


0.027 










116.38876 


-38.13771 


11.174 


0.004 


11.205 


0.005 


11.210 


0.013 


11.179 


0.018 






116.39985 


-37.77279 


12.182 


0.005 


12.145 


0.007 


12.004 


0.021 


12.057 


0.033 
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TABLE 5 — Continued 



RA 


DEC 


[3.6] 


a([3.6]) 


[4.5] 


<x([4.5]) 


[5.8 o-([5.8]) 


[8.0] 


<x([8.0]) 


[24] 


a([24]) 


J2000 


J2000 


mag 


mag 
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mag 


mag 


116.43397 


-37.73853 


12.945 


0.008 


12.946 


0.011 


12.915 


0.037 


12.959 


0.064 






116.44578 


-37.93944 


12.669 


0.007 


12.620 


0.008 


12.502 


0.036 


12.574 


0.048 






116.47578 


-37.86585 


12.499 


0.006 


12.224 


0.009 


12.217 


0.027 


12.144 


0.041 






116.48650 


-37.80289 


13.630 


0.011 


13.542 


0.013 


13.509 


0.059 


13.538 


0.089 






116.53485 


-37.89772 


11.425 


0.004 


11.422 


0.005 


11.362 


0.014 


11.341 


0.018 






116.54976 


-38.10431 


12.338 


0.006 


12.214 


0.007 


12.181 


0.025 


12.254 


0.036 






116.56485 


-38.46915 






















116.60356 


-38.01960 


10.537 


0.003 


10.612 


0.004 


10.502 


0.009 


10.481 


0.011 


10.510 


0.073 


116.61087 


-38.08359 


11.453 


0.004 


11.471 


0.005 


11.436 


0.015 


11.413 


0.015 






116.62418 


-37.87896 


12.131 


0.006 


12.064 


0.007 


12.075 


0.021 


12.005 
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116.70172 


-38.14828 


8.436 


0.002 


8.416 


0.002 


8.415 


0.003 


8.409 


0.003 


8.260 


0.031 


116.74747 


-38.37820 (#7) 


















10.154 


0.060 


116.74952 


-37.83876 


11.890 


0.005 


11.850 


0.006 


11.768 


0.018 


11.779 


0.034 






116.76556 


-37.95779 


11.713 


0.004 


11.662 


0.005 


11.580 


0.016 


11.589 


0.021 






116.76661 


-38.21130 


11.541 


0.004 


11.534 


0.005 


11.492 


0.015 


11.472 


0.031 






116.77044 


-38.22134 


11.669 


0.004 


11.734 


0.006 


11.629 


0.016 


11.649 


0.025 






116.80532 


-37.93312 


11.465 


0.004 


11.479 


0.005 


11.417 


0.014 


11.378 


0.021 






116.87226 


-37.88074 






















116.93404 


-37.98419 






















Members from Platais et al. (2001) 


115.92881 


-38.20190 


6.835 


0.002 


6.856 


0.002 


6.892 


0.002 


6.811 


0.001 






115.97642 


-38.00076 


10.367 


0.003 


10.306 


0.004 


10.274 


0.008 


10.261 


0.010 


10.113 


0.052 


116.14244 


-37.94292 


6.267 


0.002 


6.271 


0.002 


6.329 


0.001 


6.275 


0.001 


6.305 


0.031 


116.18701 


-38.05376 (#12) 


7.694 


0.002 


7.684 


0.002 


7.633 


0.002 


7.582 


0.002 


6.033 


0.031 


116.24871 


-38.09417 


7.880 


0.002 


7.916 


0.002 


7.902 


0.003 


7.856 


0.002 


7.844 


0.032 


116.54395 


-37.93366 


6.181 


0.002 


6.149 


0.002 


6.174 


0.001 


6.158 


0.001 


6.206 


0.031 


116.75196 


-37.84372 


9.880 


0.002 


9.901 


0.003 


9.851 


0.006 


9.815 


0.006 


9.814 


0.045 


Members from Hiinsch et al. (2003") 


116.02866 


-37.98943 


11.336 


0.004 


11.271 


0.005 


11.192 


0.013 


11.207 


0.018 






116.10025 


-38.59977 (#11) 


















8.961 


0.038 


116.53718 


-37.39910 (#10) 


















9.171 


0.038 



numbers next to the eoordinates designate 24 /^ni excees sourees from Table [7] 



TABLE 6 

Spitzer 3.6^M - 24^M PHOTOMETRY OF MEBERS OF NGC2451 B 



RA 
J2000 



DEC 
J2000 



[3.6] o-([3.6]) [4.5] o-([4.5]) [5.8 o-([5.8]) [8.0] cr([8.0]) [24] o-([24]) 



mag 



mag 



mag 



mag 



mag 



mag 



mag 



115.72778 
115.75046 
115.75243 
115.77086 
115.77200 
115.79673 
115.79724 
115.80643 
115.82157 



-37.95026 
-37.79300 
-37.95836 
-37.81706 
-37.98902 
-38.09564 
-38.15298 
-37.79030 
-38.28752 



TABLE 6 — Continued 



RA 


DEC 


[3.6] 


a([3.6]) 


[4.5] 


^(['1.5]) 


[5.8 


t([.-).8]) 


[8.0] 


<t([8.0]) 


[24] 


a([24]) 


.12000 


.12000 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


115.82535 


-37.76639 






















115.83157 


-38.09364 






















115.85167 


-37.93698 






















115.86620 


-37.92305 






















115.86779 


-37.94924 






















115.86904 


-38.05129 






















115.86959 


-37.94765 






















115.87277 


-37.75060 


7.541 


0.002 


7.547 


0.002 


7.515 


0.002 


7.563 


0.002 






115.87687 


-38.35557 






















115.90090 


-37.87689 


13.309 


0.009 


13.271 


0.012 


13.328 


0.047 


13.257 


0.089 






115.91875 


-38.14002 


11.166 


0.003 


11.186 


0.004 


11.152 


0.013 


11.213 


0.016 






115.93221 


-37.69007 






11.534 


0.005 






11.488 


0.020 






115.94174 


-37.91011 


13.250 


0.010 


13.170 


0.012 


13.205 


0.047 


13.166 


0.054 






115.95562 


-37.89574 (#14) 


9..343 


0.002 


9.373 


0.003 


9.310 


0.005 


9.308 


0.004 


8.549 


0.034 


115.95767 


-38.21386 


12.271 


0.005 


12.329 


0.007 


12.285 


0.025 


12.284 


0.039 






115.96495 


-38.22524 


13.292 


0.009 


13.379 


0.013 


13.277 


0.048 


13.476 


0.086 






115.98142 


-38.18655 


11.225 


0.004 


11.204 


0.005 


11.237 


0.013 


11.165 


0.016 






115.98473 


-37.77708 


12.048 


0.006 


12.052 


0.007 


11.980 


0.021 


11.922 


0.033 






115.99301 


-37.62355 






10.193 


0.003 






10.141 


0.009 


10.300 


0.048 


116.01189 


-38.12175 


12.425 


0.007 


12.425 


0.009 


12.372 


0.025 


12.417 


0.043 






116.01443 


-38.12092 


13.103 


0.008 


13.044 


0.014 


12.937 


0.032 


13.123 


0.073 






116.03135 


-37.74863 


13.200 


0.009 


13.131 


0.013 


12.974 


0.040 


13.169 


0.088 






116.04131 


-37.91467 


11.110 


0.004 


11.071 


0.004 


11.042 


0.011 


11.066 


0.015 






116.06585 


-37.76541 


11.889 


0.005 


11.840 


0.006 


11.798 


0.016 


11.778 


0.032 






116.08724 


-37.88116 


12.359 


0.006 


12.348 


0.007 


12.303 


0.025 


12.307 


0.039 






116.09396 


-37.88209 


11.148 


0.004 


11.183 


0.005 


11.116 


0.012 


11.132 


0.018 






116.09680 


-37.92455 


10.259 


0.003 


10.334 


0.004 


10.236 


0.008 


10.210 


0.008 


10.229 


0.062 


116.12193 


-37.77896 


10.950 


0.003 


10.973 


0.004 


10.929 


0.011 


10.903 


0.014 






116.13069 


-37.72335 


13.300 


0.009 


13.337 


0.012 


13.171 


0.044 


13.421 


0.081 






116.13331 


-37.74069 


13.008 


0.025 


12.923 


0.034 


12.936 


0.078 


12.797 


0.065 






116.15142 


-37.81565 


12.066 


0.005 


12.044 


0.007 


12.059 


0.024 


11.992 


0.036 






116.15163 


-37.90051 


11.229 


0.004 


11.224 


0.005 


11.138 


0.014 


11.184 


0.014 






116.15218 


-38.40723 


11.502 


0.004 






11.435 


0.016 










116.15345 


-37.91661 


12.568 


0.007 


12.576 


0.009 


12.383 


0.024 


12.582 


0.047 






116.15422 


-37.82180 


12.2.38 


0.006 


12.215 


0.007 


12.116 


0.023 


12.206 


0.046 






116.16007 


-37.86702 


14.124 


0.015 


14.029 


0.019 


14.097 


0.088 










116.17784 


-38.18400 


11.244 


0.004 


11.278 


0.005 


11.233 


0.013 


11.279 


0.016 






116.17909 


-37.68557 


13.816 


0.012 


13.784 


0.015 


13.806 


0.076 


13.748 


0.119 






116.18128 


-38.30692 


12.564 


0.006 






12.451 


0.028 










116.18606 


-37.80272 


12.694 


0.007 


12.653 


0.008 


12.619 


0.031 


12.696 


0.046 






116.20740 


-38.13024 


9.980 


0.002 


9.911 


0.003 


9.925 


0.007 


9.908 


0.007 


9.932 


0.042 


116.20831 


-38.33110 


12.683 


0.007 






12.591 


0.034 










116.21406 


-38.13371 


11.306 


0.004 


11.252 


0.005 


11.251 


0.014 


11.249 


0.016 






116.21574 


-38.16710 


12.728 


0.007 


12.709 


0.009 


12.663 


0.032 


12.553 


0.048 






116.21670 


-38.28414 


11.407 


0.004 


11.406 


0.005 


11.324 


0.013 


11.374 


0.017 






116.21722 


-38.06004 


11.936 


0.005 


11.940 


0.006 


11.905 


0.022 


11.927 


0.049 






116.21815 


-37.99622 


12.130 


0.008 


12.150 


0.008 


12.090 


0.034 


12.054 


0.050 






116.21956 


-37.96960 


13.577 


0.027 


13.396 


0.017 






13.642 


0.130 






116.23083 


-37.78915 


13.732 


0.011 


13.670 


0.017 


13.632 


0.084 


13.783 


0.126 






116.24108 


-38.02993 


13.476 


0.016 


13.386 


0.023 


13.265 


0.051 










116.26111 


-38.03118 


10.668 


0.004 


10.707 


0.004 


10.687 


0.014 


10.628 


0.038 






116.26560 


-38.10297 


10.927 


0.003 


10.948 


0.004 


10.889 


0.011 


10.936 


0.013 






116.26915 


-37.94337 


11.415 


0.006 


11.381 


0.032 


11.323 


0.031 


11.344 


0.040 






116.27896 


-38.20475 


12.141 


0.005 


12.104 


0.007 


12.039 


0.021 


12.041 


0.036 






116.28159 


-37.66232 


11.530 


0.004 


11.485 


0.005 


11.555 


0.017 


11.457 


0.026 







TABLE 6 — Continued 



RA 


DEC 


[3.6] 


a([3.6]) 


[4.5] 




[5.8 


t([.-,.8]) 


[8.0] 


<t([8.0]) 


[24] 


<t([24]) 


.12000 


.12000 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


iiiag 


mag 


116.28564 


-38.28589 


12.848 


0.007 






12.711 


0.032 










116.29492 


-38.32831 


13.914 


0.014 






13.883 


0.087 










116.30032 


-37.75469 


11.290 


0.004 


11.319 


0.005 


11.398 


0.014 


11.271 


0.018 






116.30482 


-37.83787 


10.996 


0.003 


10.998 


0.004 


11.039 


0.012 


10.957 


0.016 






116.30806 


-37.78869 


12.918 


0.008 


12.807 


0.009 


12.701 


0.032 


12.839 


0.045 






116.32416 


-37.81474 


11.223 


0.004 


11.257 


0.005 


11.217 


0.013 


11.241 


0.019 






116.32484 


-37.82308 


11.778 


0.005 


11.802 


0.006 


11.777 


0.017 


11.725 


0.024 






116.33729 


-37.83136 


10.373 


0.003 


10.353 


0.003 


10.383 


0.008 


10.319 


0.008 


10.624 


0.066 


116.34253 


-38.16618 


12.272 


0.006 


12.273 


0.007 


12.169 


0.023 


12.144 


0.035 






116.34356 


-38.06953 


13.550 


0.010 


13.512 


0.014 


13.515 


0.081 


13.531 


0.100 






116.36942 


-38.19338 


13.407 


0.009 


13.342 


0.012 


13.361 


0.057 


13.269 


0.077 






116.37029 


-38.09812 


12.924 


0.008 


12.995 


0.011 


13.111 


0.064 


12.889 


0.049 






116.38677 


-37.85640 


12.246 


0.005 


12.193 


0.007 


12.095 


0.019 


12.194 


0.048 






116.39377 


-37.86470 


11.217 


0.004 


11.271 


0.005 


11.270 


0.014 


11.227 


0.014 






116.39548 


-37.83491 


13.450 


0.012 


13.252 


0.012 


13.283 


0.056 


13.254 


0.086 






116.39838 


-37.91670 


12.237 


0.007 


12.423 


0.012 






12.508 


0.093 






116.40154 


-37.77890 


11.915 


0.005 


11.944 


0.006 


11.837 


0.018 


11.981 


0.029 






116.40282 


-37.94396 


13.175 


0.008 


13.212 


0.013 


13.662 


0.077 


13.289 


0.082 






116.40669 


-37.50750 






















116.44107 


-37.89720 


10.869 


0.003 


10.849 


0.004 


10.871 


0.011 


10.733 


0.012 


10.368 


0.049 


116.45256 


-37.54667 






9.088 


0.003 






9.057 


0.004 


8.977 


0.038 


116.46779 


-38.52593 






















116.48573 


-37.95175 


13.559 


0.010 


13.329 


0.013 


13.270 


0.054 


13.560 


0.088 






116.49277 


-37.99629 


10.329 


0.003 


10.338 


0.003 


10.336 


0.008 


10.311 


0.009 


10.326 


0.062 


116.50293 


-37.67097 


10.801 


0.003 


10.819 


0.004 


10.805 


0.010 


10.725 


0.013 






116.50903 


-37.86043 


10.913 


0.003 


10.987 


0.004 


10.915 


0.011 


10.879 


0.014 






116.50935 


-38.07950 


12.816 


0.007 


12.745 


0.009 


12.697 


0.038 


12.680 


0.050 






116.51448 


-38.01421 


13.498 


0.009 


13.419 


0.013 


13.381 


0.045 


13.579 


0.080 






116.54938 


-38.21002 


11.339 


0.005 


11.312 


0.005 


11.275 


0.014 


11.218 


0.016 


11.221 


0.127 


116.56459 


-37.74740 


12.325 


0.006 


12.312 


0.007 


12.339 


0.028 


12.239 


0.037 






116.57994 


-37.56290 






11.977 


0.006 






11.911 


0.031 






116.58564 


-37.96366 


13.892 


0.013 


13.785 


0.019 


13.788 


0.088 


13.711 


0.122 






116.60560 


-38.23373 


14.052 


0.013 


13.853 


0.016 


13.681 


0.076 


13.519 


0.083 






116.61145 


-38.10934 


12.265 


0.005 


12.198 


0.007 


12.130 


0.025 


12.139 


0.037 






116.66325 


-38.32594 


13.385 


0.009 






13.300 


0.058 










116.66850 


-38.04359 


12.217 


0.006 


12.233 


0.008 


12.193 


0.027 


12.147 


0.050 






116.67281 


-37.84615 


12.832 


0.007 


12.574 


0.008 


12.434 


0.029 


12.514 


0.043 






116.67488 


-37.98390 (#13) 


12.227 


0.005 


12.230 


0.007 


12.081 


0.021 


12.105 


0.041 


10.955 


0.072 


116.68185 


-38.28975 


11.682 


0.004 






11.726 


0.018 










116.68972 


-37.82607 


11.071 


0.004 


11.073 


0.004 


11.063 


0.011 


11.020 


0.014 






116.69822 


-38.00890 


13.726 


0.011 


13.540 


0.014 


13.406 


0.057 


13.888 


0.158 






116.72084 


-37.71950 


11.653 


0.004 


11.633 


0.006 


11.616 


0.019 


11.574 


0.022 






116.72728 


-37.59192 






11.864 


0.006 






11.802 


0.025 






116.73331 


-38.35595 


9.071 


0.003 














9.068 


0.039 


116.74249 


-37.63916 (#18) 


10.586 


0.003 


10.587 


0.004 


10.589 


0.009 


10.578 


0.010 


10.302 


0.059 


116.75339 


-37.96264 


13.379 


0.010 


13.279 


0.012 


13.266 


0.053 


13.330 


0.081 






116.76585 


-37.96487 


12.370 


0.006 


12.367 


0.008 


12.292 


0.027 


12.268 


0.043 






1 io. * / Don 


-0 ^ . t ooou 


10 111 


u.uuo 




u.uuo 




n 099 


19 1 9n 


U.Uf4U 






116.79779 


-37.79682 






















116.81857 


-38.30363 


10.169 


0.002 






10.152 


0.007 










116.84155 


-38.30055 


12.075 


0.005 






12.043 


0.022 










116.84765 


-37.71917 






















116.86294 


-38.13213 






















116.92513 


-38.18621 






















116.94233 


-38.15534 























TABLE 6 — Continued 



RA 


DEC 


[3.6] 


a([3.6]) 


[4.5] 


<x([4.5]) 


[5.8 


^([5.8]) 


[8.0] 


<x([8.0]) 


[24] 


a([24]) 


J2000 


J2000 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


116.94269 


-37.51560 






















116.94462 


-37.88126 






















Members from Hiinsch et al. (2003) 


115.95710 


-37.96688 


9.491 


0.002 


9.460 


0.003 


9.460 


0.005 


9.456 


0.006 


9.345 


0.036 


116.03583 


-37.98349 


9.044 


0.003 


9.059 


0.003 


9.044 


0.004 


9.011 


0.004 


9.073 


0.032 


116.07976 


-38.00541 


11.294 


0.004 


11.269 


0.005 


11.215 


0.013 


11.226 


0.018 






116.15953 


-38.22106 


9.010 


0.002 


8.985 


0.003 


8.932 


0.004 


8.922 


0.004 


8.839 


0.033 


116.20300 


-38.09869 


11.365 


0.004 


11.459 


0.005 


11.417 


0.015 


11.325 


0.018 






116.21267 


-37.66869 (#15) 






9.797 


0.003 


9.883 


0.007 


9.755 


0.005 


9.285 


0.034 


116.36051 


-38.19643 (#19) 


10.036 


0.003 


10.056 


0.003 


10.017 


0.007 


10.029 


0.007 


9.704 


0.038 


116.38235 


-37.50727 


















10.648 


0.069 


116.39170 


-37.79194 


11.008 


0.004 


11.017 


0.004 


10.964 


0.011 


10.953 


0.014 


10.912 


0.069 


116.50776 


-38.00752 (#16) 


10.562 


0.004 


10.540 


0.004 


10.495 


0.009 


10.474 


0.010 


9.983 


0.043 



numbers next to the coordinates designate 24 /j.ni excees sources from Table [71 
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